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Electroweak corrections to three-jet production in e + e annihilation 
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We compute the electroweak 0(a 3 a s ) corrections to three-jet production and related event-shape 
observables at electron-positron colliders. We properly account for the experimental photon isolation 
criteria and for the corrections to the total hadronic cross section. Corrections to the three-jet rate 
and to normalised event-shape distributions turn out to be at the few-per-cent level. 
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Precision QCD studies at electron-positron colliders 
rely on the measurement of the three-jet production cross 
section and related event-shape observables. The devia- 
tion from simple two-jet configurations is proportional 
to the strong coupling constant a s , so that by compar- 
ing the measured three-jet rate and related event shapes 
(see, e.g., Ref. with the theoretical predictions, one 
can determine a s . Including electroweak coupling fac- 
tors, the leading-order (LO) contribution to this process 
is of order a 2 a s . 

Owing to recent calculational progress, the QCD pre- 
dictions for event shapes @, H[ and three-jet produc- 
tion 0, [E[ are accurate to next-to-next-to-leading order 
(NNLO, a 2 al) in QCD perturbation theory. Depend- 
ing on the observable under consideration, the numeri- 
cal magnitude of the NNLO corrections varies between 
three and twenty per cent. Inclusion of these corrections 
results in an estimated residual uncertainty of the QCD 
prediction from missing higher orders at the level of be- 
low five per cent for the event-shape distributions, and 
below one per cent for the three-jet cross section. At this 
level of theoretical precision, higher-order electroweak ef- 
fects could be of comparable magnitude. At present, only 
partial calculations of electroweak corrections to three- 
jet production and event shapes are available [f|, which 
can not be compared with experimental data directly. In 
this work, we present the first calculation of the NLO 
electroweak (a 3 a s ) corrections to three-jet observables 
in e + e~ collisions including the quark-antiquark-photon 
(qq~f) final states. Note that the QCD corrections to 
these final states are of the same perturbative order as 
the genuine electroweak corrections to quark-antiquark- 
gluon (qqg) final states. Since photons produced in as- 
sociation with hadrons can never be fully isolated, both 
types of corrections have to be taken into account. 

Event-shape measurements at LEP usually rely on a 
standard set of six variables y, defined for example in 
Ref. thrust T, C-parameter, heavy jet mass p, wide 
and total jet broadenings By/ and B^, and two-to-three- 
jet transition parameter in the Durham algorithm I3. 
The experimentally measured event-shape distribution 

1 da 



is normalised to the total hadronic cross section. In the 
perturbative expansion, it turns out to be most appropri- 
ate to consider the expansion of this ratio, which reads 
to NNLO in QCD and NLO in the electroweak theory 
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where the fact is used that the perturbative expansion 
of Chad starts at order a 2 . The calculation of the QCD 
coefficients A, B, and C is described in Refs. [2|, |3(. The 
LO purely electromagnetic contribution <5 7 arises from 
tree-level qqj final states without a gluon. The NLO 
electroweak coefficient Sew receives contributions from 
the 0(a) correction to the hadronic cross section, 



Chad = Cq 
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and from the genuine <D(a 3 a s ) contribution to the event- 
shape distribution 
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yields the full NLO electroweak correction. Both terms 
are to be evaluated with the same event-selection cuts. 
As shown in the following, many of the numerically domi- 
nant contributions, especially from initial-state radiation, 
cancel in this difference. 

In the experimental measurement of three-jet observ- 
ables at electron-positron centre-of-mass energy yfs, sev- 
eral cuts are applied to reduce the contributions from 
photonic radiation. In our calculation, we apply the 
criteria used in the ALEPH analysis Q. Very similar 
criteria were also applied by the other experiments [![. 
Particles contribute to the final state only if they are 
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within the detector acceptance, defined by the produc- 
tion angle relative to the beam direction, |cos 6>| < 0.965. 
Events are accepted if the reconstructed invariant mass 
squared s' of the final-state particles is larger than s cut = 
0.81s. To reduce the contribution from hard photon ra- 
diation, the final-state particles are clustered into jets 
using the Durham algorithm with resolution parameter 
Ucut = 0.002. If one of the resulting jets contains a pho- 
ton carrying a fraction z 7 > z 7tCU t = 0.9 of the jet energy, 
it is considered to be an isolated photon, and the event is 
discarded. The event-shape variables arc then computed 
in the centre-of-mass frame of the final-state momenta, 
which can be boosted relative to the e + e~ centre-of-mass 
frame, if particles are outside the detector acceptance. 

In the computation of the 0{a) corrections to the to- 
tal hadronic cross section, we include the virtual elec- 
troweak corrections to qq final states, and the real radia- 
tion corrections from qqj final states, provided the above 
event-selection criteria are fulfilled. The corrections to 
the event-shape distributions receive contributions from 
the virtual electroweak corrections to the qqg final state, 
the virtual QCD corrections to the qqj final state, and 
from the real radiation qqg^ final state. To separate the 
divergent real radiation contributions, we used both the 
dipole subtraction method d, [l(| and phase-space slic- 
ing [Til ] . resulting in two independent implementations. 
Soft singularities arc present in the virtual and real cor- 
rections. They are regularized dimensionally or with in- 
finitesimal photon and gluon masses, and cancel in the 
sum. Collinear singularities from photon radiation off 
the incoming leptons (initial-state radiation, ISR) are 
only partially cancelled. The left-over collinear ISR sin- 
gularity is regularized by the electron mass and absorbed 
into the initial-state radiator function, which we consider 
either at fixed order, or in a leading-logarithmic (LL) 
resummation [T3 ]. Owing to the specific nature of the 
event selection, also photon radiation off the outgoing 
quarks (final-state radiation, FSR) is only partially can- 
celled. The left-over FSR singularity arises from the iso- 
lated photon definition, which vetoes on photon jets with 
z-y > z 7iCUt . This singularity is absorbed into the photon 
fragmentation function, which we apply in the fixed-order 
approach of Ref. [l3| . For the non-perturbative contribu- 
tion to this function, we use the 0(a) two-parameter fit 
of ALEPH 0. The fragmentation contribution derived 
in Ref. [l3| is based on phase space slicing and dimen- 
sional regularization. We recomputed this contribution 
using subtraction and mass regularization [Toj . 

The Feynman diagrams for the virtual corrections are 
generated with FeynArts [HI, Ojj]. Using two inde- 
pendent inhouse Mathematica routines, one of which 
builds upon FormCalc jl7| . each diagram is expressed 
in terms of standard matrix elements and coefficients of 
tensor integrals. The tensor integral coefficients are nu- 
merically reduced to standard scalar integrals using the 
methods described in Refs. [l|| [2(|. The scalar mas- 
ter integrals are evaluated using the methods and results 
of Refs. [U [22|, HH, where UV divergences are regular- 



ized dimensionally. For IR divergences two alternative 
regularizations are employed, one that is fully based on 
dimensional regularization with massless light fermions, 
gluons, and photons, and another that is based on in- 
finitesimal photon and gluon masses and small fermion 
masses. The loop integrals are translated from one 
scheme to the other as described in Ref. 24]. 

The Z-boson resonance is described in the complex- 
mass scheme [2^, [26| , and its mass is fixed from the com- 
plex pole. The electromagnetic couplings appearing in 
LO are parametrized in the scheme, i.e., they are 
fixed via 

a = a Gfi = V2G^M^ (l - M&/Aff ) /tt. 

As the leading electromagnetic corrections are related 
to the emission of real photons, we fix the electromag- 
netic coupling appearing in the relative corrections by 
a = a(0), which is the appropriate choice for the lead- 
ing photonic corrections. Accordingly the cross section 
for e + e _ — * qqg is proportional to Qg a s while the elec- 
troweak corrections to this process are proportional to 
a(Q)o? Gij a s . 

We performed two independent calculations of all in- 
gredients resulting in two independent Fortran codes, 
one of them being an extension of Pole 18|. 

We use the following values of the electroweak and 
QCD parameters: 

G M = 1.16637 x 10~ 5 GcV~ 2 , M H = 120 GeV 
a(0) = 1/137.03599911, a s (M z ) = 0.1176, 

to c = 0.51099892 MeV, m t = 171.0 GeV. (4) 

Because we employ a fixed width in the resonant W- and 
Z-boson propagators in contrast to the approach used at 
LEP to fit the W and Z resonances, where running widths 
are taken, we have to convert the "on-shell" values of 
M£ EP and P^ EP (V = W, Z), resulting from LEP, to the 
"pole values" denoted by My and Ty, leading to [27| : 

A% = 80.375 . . . GeV, T w = 2.140 . . . GeV, 
M z = 91.1535... GeV, T z = 2.4943 ... GeV. (5) 

In the final state we take all light quarks into account, 
including b quarks. 

In Fig. [TJ we display the total hadronic cross section 
Chad for the above event-selection criteria including NLO 
electroweak corrections and the relative corrections sep- 
arately. For the latter, "full" and "weak" refers to the 
electroweak NLO corrections with and without purely 
photonic corrections, respectively, and "h.o. LL" indi- 
cates the inclusion of the higher-order ISR effects. For 
most energies, the full 0(a) corrections are sizable and 
negative, ranging between —30% at the Z peak and about 
— 10% at energies above and below. The numerically 
largest contribution is always due to ISR. However, above 
the Z resonance up to about 110 GeV the corrections are 
positive and of the order of the Born cross section due 
to the well-known radiative return phenomenon [28j that 
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FIG. 1: Total hadronic cross section Chad 
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FIG. 2: Differential thrust distribution at y/s = Mz- 
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FIG. 3: Differential thrust distribution at ^fs — Mz nor- 
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FIG. 4: Three-jet rate at y/s = Mz normalised to (Thad- 



occurs in this region because of our choice of s cu t- As it is 
not relevant experimentally, we did not fully resolve it in 
Fig. [T] Below 60 GeV and above 120 GeV the magnitude 
of the corrections is increased due to LL rcsummation 
of ISR, whereas it is decreased in the region in between. 
The virtual one-loop weak corrections (from fermionic 
and massive bosonic loops) yield only a moderate cor- 
rection between —5 and +5%. The increase for energies 
above 1 TeV is due to electroweak Sudakov logarithms. 

Using the same event-selection cuts, Fig. [5] displays the 
differential thrust distribution at y/s = Mz, including 
NLO electroweak contributions. The distributions are 
weighted by (1 — T), evaluated at each bin centre. The 
Born contribution is given by the A-term of ([1]), while the 
full 0(a) corrections contain the tree-level qqj contribu- 
tion (5 7 and the NLO electroweak contribution &g. Again, 
we observe large negative corrections due to ISR, and 
moderate weak corrections. The corrections are largely 
constant for T < 0.95, where the isolated photon veto re- 
jects all contributions from qqj final states. For T > 0.95, 



corresponding to the two-jet limit, we find a substantial 
contribution from qq^ final states already at LO (a 3 ). 
Moreover, it turns out that the electromagnetic correc- 
tions depend non-trivially on the event-selection cuts. 

In expanding the corrections according to (fTJ), and re- 
taining only terms up to LO in a s , we obtain the gen- 
uine electroweak corrections to normalised event-shape 
distributions, which we display for thrust at i/s = Mz 
in Fig. [3] Again, the Born contribution is given by the 
A-term, while the 0(a) corrections now consist of <5 7 and 
(Jew- It can be seen very clearly that the large ISR cor- 
rections cancel between the event-shape distribution and 
the normalisation to the hadronic cross section, resulting 
in electroweak corrections of a few per cent. Moreover, 
effects from ISR resummation are largely reduced as well. 
The purely weak corrections are below 0.5 per mille. 

As shown in Fig. 01 we observe a similar behaviour 
for the three-jet rate in the region y cut > 0.002. For 
Vcut ^$ 0.002, qq^f final states contribute and lead to a 
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distortion of the shape. For y cut < 0.0005 the three-jet 
rate becomes larger than Chad- I n this region, fixed order 
perturbation theory becomes unreliable due to large log- 
arithmic corrections at all orders. By improving the LO 
QCD prediction used here with NLO and NNLO QCD 
corrections, which are large and negative [H, 0] for small 
2/cut) it is possible to extend the range of validity of the 
fixed order predictions. 

In Ref . @ , another calculation of electroweak corrections 
to three-jet observables was performed, which differs in 
two important aspects from the work presented here. It 
considered only the corrections to qqg final states, while 
qq'y final states at LO and NLO were not taken into ac- 
count. To remove singularities associated with infrared 
gluons in qq^/g final states, event-shape observables were 
calculated from the reconstructed jet momenta and not 
from the parton momenta, as used in experiment and 
in our work. Moreover, the NLO electroweak correc- 
tions to the hadronic cross section were not taken into 
account, such that only unnormalised distributions were 
considered. Owing to these substantial differences, a di- 
rect comparison with the results of Ref. Q is not possi- 
ble. Taking care of the different renormalizaion of a, we 



do observe, however, in the unnormalized distribution, 
Fig. [21 that the relative size of the O(a) weak and exact 
corrections, and of the LL-improved corrections to the 
thrust distribution agree at the per-cent level with the 
results of Ref. [||, except in the region (1 — T) < 0.05, 
where qq^ final states contribute. 

Data on event-shape distributions and jet cross sec- 
tions have been corrected for photonic radiation effects 
modelled by standard LL parton-shower Monte Carlo 
programs. They can thus not be compared directly with 
the NLO electroweak corrections computed here. Incor- 
poration of these corrections requires a more profound 
reanalysis of LEP data, in order to quantify the impact 
of the NLO electroweak corrections on precision QCD 
studies, such as the precise extraction of the strong cou- 
pling constant at NNLO in QCD [29| . 
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